Existence and uniqueness for a class of nonlinear 
higher-order partial differential equations in the complex 

plane 

O. Costin* and S. Tanveer 
February 1, 2008 



Abstract 

We prove existence and uniqueness results for nonlinear third order partial differential 
equations of the form 

3 

!t-fy yy = Y, h Ay, t ;f) f U) +r(y,t) 

0=0 

where superscript j denotes the j-th partial derivative with respect to y. The inhomogeneous 
term r, the coefficients bj and the initial condition f(y, 0) are required to vanish algebraically 
for large \y\ in a wide enough sector in the complex y- plane. Using methods related to Borel 
summation, a unique solution is shown to exist that is analytic in y for all large \y\ in a 
sector. Three partial differential equations arising in the context of Hele-Shaw fingering 
and dendritic crystal growth are shown to be of this form after appropriate transformation, 
and then precise results are obtained for them. The implications of the rigorous analysis on 
some similarity solutions, formerly hypothesized in two of these examples, are examined. 



1 Introduction 

The theory of partial differential equations (PDEs), when one or more of the independent vari- 
ables are in the complex plane, appears to be largely undeveloped. The classic Cauchy-Kowaleski 
theorem holds for a system of first-order equations (or those equivalent to it) when the quasi- 
linear equations have analytic coefficients and analytic initial data is specified on an analytic 
but non-characteristic curve. Then, the C-K theorem guarantees the local existence and unique- 
ness of analytic solutions. As is well known, its proof relies on the convergence of local power 
series expansions and, without the given hypotheses, the power series may have zero radius of 
convergence and the C-K method does not yield solutions. 
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Relatedly, not much is known in general for higher-order nonlinear partial differential equa- 
tions in the complex plane, even for analytic initial conditions and analytic dependence of the 
coefficients. The only work we are aware of on nonlinear partial differential equations in the 
complex plane involving higher spatial derivatives is that of Sammartino & Cafiisch pi ^ who 
proved among other results the existence of a solution to nonlinear Prandtl boundary layer 
equations for analytic initial data. This work involved inversion of the heat operator dt — dyy 
and using the abstract Cauchy-Kowalewski theorem for the resulting integral equation. Unfor- 
tunately, this methodology cannot be adapted to our problem. The coefficients of the highest 
(third order) spatial derivatives in our equation depend on the unknown function as well. These 
terms cannot be controlled by inversion of a linear operator and estimates of the kernel, as used 
by Sammartino & Cafiisch. Instead, the essence of the methodology introduced here is the use 
of large y asymptotics, conveniently expressed in terms of the behavior of the unknown function 
in the Borel transform variable p for small p. The choice of appropriate Banach spaces proves 
to be crucial, and after this choice the contraction mapping argument itself is not difficult. 

One aim of the present paper is to obtain actual solutions with good smoothness and asymp- 
totic properties for a class of PDEs, when power series solutions may have zero radius of conver- 
gence. Our approach based on Borel summation techniques provides at the same time appropri- 
ate existence and uniqueness results for a class of nonlinear PDEs in the complex domain. 

Keeping in mind applications, we develop the framework for certain higher order partial dif- 
ferential equations in a domain where one of the independent variables (y in this case) is complex, 
while the other (t) is real. While more sweeping generalizations are under way, the current paper 
is restricted in scope by the applications we have in mind and simplicity of exposition. 

There is a class of nonlinear PDEs that have recently arisen in applications. The basic features 
of the application problems is that in the absence of a regularization (like surface tension) the 
initial value problem in the real domain is relatively simple; yet ill-posed in the sense of Hadamard 
for any Sobolev norm on the real domain. However, the analytically continued equations into the 
complex spatial domain are well-posed, even without a regularization term. Earlier, Garabedian 
p] had recognized the conversion of ill-posed elliptic initial value problem into a well-posed one 
by excursion into the complex plane in the spatial variable. Moore Cafiisch & Orellana 

@ ; [ft ; Cafiisch & Semmes || and Cafiisch et al Jl(| have studied solutions to the complex plane 
equations that arose from simplifications of vortex sheet evolution (in fluid mechanical contexts). 
The initial value problem in these cases is ill-posed in the real domain, though well-posed in an 
appropriate class of analytic functions in a domain in the complex plane. Study of the complex 
equations proved useful since evidence 1 1 of finite time singularities in the real domain can be 
traced to earlier singularity formation in the complex domain. 

In the physical context of Hele-Shaw dynamics, it was suggested Jl2| that it is fruitful to study 
the complex plane equations even when the initial value problem is well-posed in the physical 
domain through the addition of a small regularization term. The advantage of this procedure is 
that one can study small regularization effects by perturbing about the relatively simpler but well- 
posed zeroth order problem. The ill-posedness of the unregularized problem in the real domain, 
shown earlier by Howison Jl5[ , is transferred into ill-posedness of the analytic continuation of 
initial data to the complex plane. However, when analytic initial data is specified in a domain 
in the complex plane such as to allow for isolated singularities, there is no ill-posedness of the 
zeroth order approximation of the dynamics. This provides the basis for a perturbative study 
that includes small but nonzero regularization effects in the real domain. Consideration of an 
ensemble of complex initial conditions, subject to appropriate constraints on its behavior on the 



2 



real axis, provide a way to understand the robust features of the dynamics when regularization 
effects are small. Indeed this procedure has yielded information about how small surface tension 
can singularly perturb a smooth solution of the un-regularized dynamics jll| . It has given scaling 
results on nonlinear dendritic processes as well |TjJ. However, much of the results derived so far 
are purely formal and rely fundamentally on the existence and uniqueness of analytic solutions to 
certain higher order nonlinear partial differential equations in a sector in the complex plane, with 
imposed far-field matching conditions. Indeed, in a more general context, one can expect that 
whenever regularization appears in the form of a small coefficient multiplying the highest spatial 
derivative, the resulting asymptotic equation in the neighborhood of initial complex singularities 
will satisfy higher order nonlinear partial differential equation with sectorial far-field matching 
condition in the complex plane of the type shown in Examples 1-3. Hence, there is a need to 
develop a general theory in this direction. 

2 Problem Statement and Main Result 

We seek to prove the existence and uniqueness of solutions f(y,t) to the initial value problem 
for a general class of quasilinear partial differential equations of the form: 

3 

/* - fyyy = E hiv^f) f U) + r(y,t) with /(„,0) = fi(y) (1) 
i=o 

where the superscript V' refers to the j-th derivative with respect to y. The inhomogeneous 
term r(y, t) is a specified analytic function in the domain 

V po = {(y,t) :arg ye (-2^/3, 2vr/3), \y\ > Po > 0, < t < T} 

and it is assumed that in T> Po there exist constants a r > 1 and A ri with only A r allowed to 
depend on T, such that 

\y a " r(y,t)\ < A r {T) (2) 

Further in (|l]), the coefficients bj may depend on the solution /- this is how nonlincarity in the 
problem arises. It is possible to extend the current theory to include dependence of bj on fU) 
as well, though for simplicity we will restrict only to dependence on /. Further, we restrict to 
the case where each bj is given by a convergent series 

oo 

bj{ y ,fj) = J2hk(y,t) f k (3) 

fc=0 

for known bj^, analytic for y in T> Po . It will be assumed that in this domain, there exists some 
choice of positive constants f3, (Xj, and Aj,, independent of j and k (with (3 and Oj independent 
of T as well), such that 

\y a ' +kfi b jtk \ < A b (T) (4) 
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Further, the series (||) converges in the domain T>$ jP , defined as 



£V, P = {(V,t) : argy G (~ - <j>, ^ + ^ , |y| > p > p 



whclV I ■ /j. (5) 



7T 

G' 



if 

l/l < / (6) 



Condition 1 The solution f(y,t) we seek for fijy is required to be analytic for complex y in 
^(p.p f or some p > (to be determined later). In the same domain, the solution and the initial 
condition fi{y) must satisfy the condition 

\y a >- f(y,t)\ < Af(T) (7) 

for some Af that can only depend on T for (y, t) € p ■ 

It is clear that for large y such a solution / will indeed satisfy (||), the condition for the conver- 
gence of the infinite series in (jl]). The general theorem proved in this paper is that 

Theorem 2 For any T > and < <fi < tt/6, there exists p such that the partial differential 
equation (^j has a unique solution f that is analytic in y and 0(y^ 1 ) as y — > oo for (y, t) G T^p,4>- 
In fact, we have for this solution f = 0(y~ ar ) as y — > oo. 

We note that uniqueness requires analyticity and decay properties of y in a large enough sector. 
The proof of theorem (^) will have to await some definitions and lemmas. It is to be noted that, 
from a formal argument if / is small, the dominant balance for large y is between f t on the left 
of (Q) and r(y,t) on the right, indicating that f(y,t) ~ fi{y) + J r {y,t) dt. Since each of 
fi(y) and r(y,t) decays algebraically as y — > oo within T> Po at a rate y~ ar , which for a r > 1 
is much less than y _1 , this suggests suggests that other terms in the differential equation ([!]) 
should not contribute. This is in fact shown rigorously for \y\ large with arg y € (— 2tt/3, 2tt/3). 
As shall be seen in the examples, this behavior for solution / is not valid outside this sector, 
where in general, one can expect infinitely many singularities with an accumulation point at oo. 



3 Inverse Laplace Transform and Equivalent Integral Equa- 
tion 

The Inverse Laplace Transform G(p, t) of a function g(y, t) analytic in 2?0 lP and vanishing alge- 
braically as \y\ — > oo is given by: 

G(p,t) = [C-'ig}}^) = jf g(y,t)dy (8) 
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where Co is a contour as in Fig. 1 (or deformations thereof), entirely within the domain T><p tP . 
We restrict p to the domain: 

S<f, = {p : arg p e (-0,0), < \p\ < oo } 

It is easily seen that if g(y,t) = y~ a for a > 0, then G(p,t) = p a ~ 1 /T(a). From the 
following Lemma, it is clear that the same kind of behavior for the Inverse Laplace Transform 
G(p, t) can be expected for small p in 5^, with a y~ a behavior of g at oo. 



Lemma 3 If g(y,t) is analytic in y in V^.p, and satisfies 



\y a g(y,t)\ < A(T) (9) 

for a > a® > 0, then for any 5 € (0,0) the Inverse Laplace Transform G = C~ 1 g exists in S^^s 
and satisfies 

\G( P ,t)\ < c^br 1 ^ (io) 

for some C — C(S, a ). 

Proof. We first consider the case when 2 > a > ao- Let C Pl be the contour Cd in Fig. 1 
that passes through the point p\ + and given by s = p\ + \p\~ x +irexp(«0signum(r)) with 

r e (—00,00). Choosing 2p > pi > (2/\/3)p, we have \s\ > p along the contour and therefore, 
with arg(p) = 6 € (—0 + 6, — 5), 

\g(s,t)\ < A{T)\s\- a and |e sp | < e" l|p|+1 e - |r||p| sin 

Thus 



'Cpi 



e sp g( Sl t)ds 



<2A(T)e pl ^ +1 / \p 1 + \p\- 1 +ire l(t> y a e-^ rs{n5 dr 
Jo 



< KA(T)e pM \ Pl + IpI' 1 ^ 01 / e-l p l rsin<5 rfr < ^- 1 |p|"- 1 e 2 " |p| (11) 

Jo 

where K and K are constants independent of any parameter. Thus, the Lemma follows for 
2 > a > ao, if we note that T(a) is bounded in this range of a, with the bound only depending 
on a - 

For a > 2, there exists an integer k > so that a — k € (1, 2]. Taking (k — l)\h(y,t) = 
J^ o g(z,t)(y — z) k ~ 1 dz (clearly h is analytic in T>^, p and h^ k \y,t) = g(y,t)), we get 



h(y,t) 



(-y) k 



(k 



-ML [ 
- 1)! Ji 



g{mt)ip-l) k - 1 dp 



k^.k — a. 



(fc-1)! 



/OO 
A(yp, t)p- 



"(p- l)* _1 dp 
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with \A(yp,t)\ < A(T), whence 

.,, A(T)T(a-k) 

\h(y,t)\ < r i L, , 
I 71 |y|«-*r(a) 

From what has been already proved, with a — k playing the role of a, 

\c- l {h}{p,t)\ < c^^br*- 1 ^ 

T{a) 

Since G(p,t) = ( — l) k p k C~ 1 {h}(p,t), by multiplying the above equation by \p\ k , the Lemma 
follows for a > 2 as well. 
□ 

Comment 1: As mentioned before, when g(y,t) — A(t)y~ a we have G(p,t) = p a ~ x . 

In this case, the exponential factor in ([h]) can be omitted because of the algebraic behavior of 
g(y,t) for all y. This result is relevant for Examples 1-3. 

Comment 2: The constant 2p in the exponential bound can be lowered to anything exceeding 
p, but ( |lo| ) suffices for our purposes. 

Comment 3: Corollary [| below implies that for any p 6 5^, the Inverse Laplace Transform 
exists for the specified functions r(y,t), bj t k(y,t), as well as the solution f(y,t) to whose 
existence is shown in the sequel. 

Comment 4: Conversely, if G(p, t) is any integrable function satisfying the exponential bound 
in ([l(]), it is clear that the Laplace Transform along a ray 

te 

CgG = / dp e- p y G(p,t) (12) 



exists and defines an analytic function of y in the half-plane 5i [e l9 y] > 2p for 9 G (—<t>, <t>)- 
Comment 5: The next corollary shows that there exist bounds for C {bj^} and £ _1 {r} 
independent of arg p in S$, because of the assumed analyticity and decay properties in the 
region T> Po , which contains 2?0 )P . 



Corollary 4 The Inverse Laplace Transform of the coefficient functions bj^k and the inhomoge- 
neous function r(y,t) satisfy the following upper bounds for any p € 6>0 



\R(p,t)\ < ^ A r (T) \p\^~ l e 2 '°W (14) 



G 



Proof. From the assumed conditions we see that bj t k is analytic in y over T>$ ltPo for any </>i 
satisfying n/6 > 4>\ > <f>. So Lemma |] can be applied for g(y,t) — bj,hi with tf>i = <j> + 
(tt/6 — </>)/2 replacing 0, and with (5 replaced by <f>\ — (f> — (tt/6 — </>)/2, and the same applies to 
t), leading to ( p"3| ) and (JT^). In the latter case, since a r > 1, ao in Lemma |3| can be chosen 
to be 1. Thus, one can choose C2 to be independent of a r , as indicated in (|l4|). □ 

The formal Inverse Laplace Transform of (fll) with respect to y is 



3 00 

F t + p 3 F = £ (-iy £ [ B j>k * (^F) * F* k ] (p,t) + R{ P ,t) (15) 

j=0 k=0 

where the symbol * stands for convolution (see also ||). On formally integrating (|l5|) with 
respect to t, we obtain the integral equation 



F M = EE / (-l) j e-^- T) [(p>'F) * B hk * F* k ] (p,r)dr + F (p,t) 

= NF{ Pl t) (16) 

where 

F (p,t) = e-P 3t Fj(p) + f e -P 3 ^ R{p,r) dr (17) 

Jo 

Here Fj =C~ 1 {fi}- 

Our strategy is to reduce the problem of existence and uniqueness of a solution of (^) to the 
problem of existence and uniqueness of a solution of ([If]), under appropriate conditions. 

4 Solution to the Integral Equation 

To establish the existence and uniqueness of solutions to the integral equation, we need to 
introduce an appropriate function class for both the solution and the coefficient functions. 

Definition 5 Denoting by 5^ the closure of S$, dS$ — S^\S^, and K, — S$ x [0,T], we define 
for v > (later to be taken appropriately large) the norm \\ ■ ||„ as 

\\G\\ U = M sup (l + b| 2 ) e-"'*' \G(p,t)\ (18) 
(p,t)eK 



where 



[ 2(1 + s 2 ) (ln(l + s 2 ) + s arctans) | , s 

Mo -sup ^ A \ 2 ' ^=3.76... 19 
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Definition 6 We now define the following class of functions: 
A& — {F : F(-, t) analytic in and continuous in 

for t e [0,T]s.t. \\F\\„ < 00} 

It is clear that A$ forms a Banach space. 

Comment 6: Note that given G 6 A$, g(y,t) = Cq{G} exists for appropriately chosen 9 
when p is large enough so that p cos(# + arg y) > v, and that \y g(y,t)\ is bounded for y on 
any fixed ray in T>^ p . 

Lemma 7 For v > 4po + a r , Fj = C~ 1 {fi} satisfies 

\\Fj\\ v < C^)A fl {v/2)-^ +1 
while R = C~ x r satisfies the relation 

\\R\\ V < C(ct>)MT)(vm- ar+1 

and therefore 

\\F \\u < C(cf>)(T A r + A fl ){v/2)-^ +1 (20) 

PROOF. First note the bounds on R in Corollary [|. We also note that a r > 1 and that for 
v > 4po + ctr we have 



SU pi^ e -("-^)W < {ar I ( ir :~ e-^(,-2p a )-^ +1 
P r(a r ) r(a r ) 

< Ka- 1 / 2 (v/2)- a r +1 

where K is independent of any parameter. The latter inequality follows by accounting for 
Stirling's formula for T(a r ) for large a r . Similarly, 



p i-(ar) f(a r ) 



Using the definition of the v— norm and the two equations above, the inequality for follows. 
Since fi{y) is required to satisfy the same bounds as r(y, t), a similar inequality holds for ||.Fj||„. 
Now, from the relation (|l7|), 

\F (p,t)\ < \Fi(p)\ + T sup \R(p,t)\ 

0<t<T 

Therefore, @) follows. □ 

Comment 7: Not all Laplace-transformable analytic functions in 2?0 iP belong to A$. For the 
applications we have in mind, the coefficients are not bounded near p = and hence do not 
belong in A^. It is then useful to introduce the following function class: 
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Definition 8 

H = {if : H(p,t) analytic in <S , \H(p,t)\ < C|p| Q " 1 e p|p| | 
for some positive constants C and a and p which may depend on H . 



Lemma 9 If H G TL and F G A^, then for v > p + 1, H * F belongs to A^, and satisfy the 
following inequality^ 

\\H*F\\<\\\H\*\F\\\ V < CT{a){v-p)- a \\F\\ V (21) 

Proof. From elementary properties of convolution, it is clear that H * F is analytic in S$ and 
is continuous on S^, Let 9 — arg p. It is to be noted that 

\p\ 

se lH )\ds 



But 
And 



\H*F(p)\ < ||£T|*|i?|(p)|< [ \H(se ie )\\F(p 

Jo 

\H(se i8 )\ < Cs Q " 1 e |s|p 



/ S ^e^\F{p-se l<) )\d S <\\F\\ v e^\pr ^ds (22) 

Jo Jo M (i + \p\ 2 {i~ s y) 

If \p\ is large, noting that v — p > 1, we obtain from Watson's lemma, 

p v \p\ 

s <*-i e WP\F(p - se lS )\ds < KT{a)\\F\\ v —— —\v - p\~ a (23) 



I 



o 



io M (l + H 2 ) 

Now, for any other |p|, we obtain from (p2|), 

s °-W\F(p ~ se l6 )\ds < K\v - p\- a \\F\\ v - / ' 

Thus the relation (|23j) must hold in general as it subsumes the above relation when \p\ is not 
large. From this relation, ( pl| ) follows from applying the definition of \\.\\ v . 
□ 



Corollary 10 For F G A<j>, and v > 4 p + 1, 

\\B j , k *F\\ v <\\\B j>k \*\F\\\ v < KCA^aj) {v /2)~ k(3 ~ a i A b (T) \\F\\ V 

PROOF. The proof follows simply by using Lemma ||, with H replaced by Bj t k and using the 
relations in Corollary 0. □ 



1 In the following equation, || ■ ||„ is extended naturally to continuous functions in K. 
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Lemma 11 For F, G € A$ and j > 0, 



\{p> F)*G(p,t)\ < -^-^L- \\F\\ V \\G\\ V 



M (l + \p\< 



Proof. Let p = \p\ e . Then 



\(p> F)*G(p,t)\ 



S j F(s) G(p - s) ds 



< / H ds \F(s e l6 )\ \G( P - se w )\ 
Jo 



Using the definition of || • ||„, the above is bounded by 



Se^l||F|U ||G|U 



M 2 



ds 



< 



o (l + s 2 )[l + (M-s) 2 ] 

e flpl 



M (l + b| 2 ) 
The last inequality follows from the fact that 



pI 



1 



/„ (i + s 2 )[i + (H-.s) 2 ] 

and the definition of Mo in (591). □ 



= 2 



11*11, l|G||„ 



In (|p| 2 + 1) + \p\ tan" 1 \p\ 



\P\ 



4) 



Corollary 12 WWi t/ie convolution *, is a Banach Algebra and furthermore, 

\\F*G\\ V < \\F\\ U \\G\\ V 

PROOF. This follows by applying Lemma |ll| for j = and using the definition of | 
□ 



Lemma 13 For v > 2/?o + 1, 



/ (p>F) * B j>k * F* k e - p3{t - T) dr 
Jo 



< 



C{4>) 



M (l + H 2 , 

where the constant C is independent of T , but depends on 



\\\Bj,k\*\F\ h \\F\\ k v e»W T^/ 3 
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PROOF. For k > 1, from Lemma |Tl|, with G = (B j}k * F) * F*^ 1 ), we obtain 



\(p> F) * B it * F' k \ < J^L_|| Biik ,F| U | J?C 



For fc = 0, we note that 



\&F)*Bi, \ < W 11^1*1^,011 < M^Tg2) lll^.o|*mi|, 



Therefore, ([28|) holds for k = as well. Thus, for any k > 0, 



/ (p*F) * Bj, k * F* k e - p3{t - T) dr 
Jo 



\\\B jtk \*\F\\U \\F\\t f e 
Jo 



<-M^—^- \\\B,,\*\F\\\„ \\F\\ k . I e -IH 3 cos(3e)( t - 



M (l+W 2 ) 

where p = |p|e ie . On integrating with respect to r, we obtain 



□ 



Definition 14 For H e H, F and h in A$, define ho — and for k > 1, 

h k = H * [(F + h)* k — F* k }. 

Lemma 15 For v > p+1. 

k-i 



\\h k \\u < k (||F||„+||fc||„) \\\H\*\h\\\ u 



Proof. We prove this by induction, k — follows trivially since ho = 0. The case of k 
obvious from f|29|). Assume the formula (|30|) holds for all k < I. Then 



\\hi+i\\ v = \\H * (F + h) * (F + h)*' l -H*F*F* l \\ v = \\H * h * (F + h)* 1 + F * h t 
On using ( |30| ) for k = I, we get 

< iii^i * hiu mu + + in(\\F\u + * nil 

<(! + l)(\\F\\„ + \\h\\ v ) 1 \\\H\*\h\\\ v 
Thus the formula (Sol) holds for k = I + 1. □ 
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Lemma 16 For F and h in and v > 2po + 1, and k > 1, 

(p> [F + h]) * Bj t k * (f + /i) * - (p>F) * B jtk * F* 

(\\F\\ v + \\h\\«) {*||^|U|||B Jlfc | 



< 



Mo(l + |p| 2 ) 



\\B jtk *(F + h)\\ v \\h\\ v } (31) 



Proof. It is clear that 



{p> [F + h]) *B Jik *(F + h) - (p>F) * B hk * F 



< 



{p J h)*B 3 , k * F + /i 



+ \(p>F)*h k \ (32) 



where H is now replaced by Bj tk in the definition of h k in (|29|). Applying Lemma [II] and 
Corollary \L2 to the first term, we obtain for k > 1 



ip*h) * B jik * (f + hj 



< 



\p\i e "\p\ 



■\\h\M\B jtk *(F + h)\\J\\F\\ v + \\h\ 



k — l 



M (l + H 2 ) 

On the other hand, applying Lemma [ll] and Lemma |lR with H = Bj tk and p replaced by 
2po, we obtain 

\{p>F)*h k \ < ^ + k\\F\\„\\B jtk * h\\ u [\\F\\ v + 

Combining the previous two equations, and using it in (|3^), we obtain the proof of Lemma fL6| 
by noting that \\B j;k * h\\ v < \\\B jtk \ * \h\\\ u . □ 



Lemma 17 For v > 2po + 1, 



{(^[F + h}) * B jfi - {p>F) * B jt0 } e -f 3 (*-) dr 



< 



M (l + H 2 ) 



lll%o|*NIL (33) 



Proof. We note that 

| {jP [F + h])* Bjfi - ( P 3 F) *B jt o\ < | { P J h) * B jfi | 



< 



A/o(l + bl 2 ) 



2Tll|Bi,o|*Nlk 
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Further, as before in the proof of Lemma O 



-V\t-r) 



Combining the two equations above, the Lemma follows. 



□ 



Lemma 18 For v > 2po + 1, and k > 1, 



{(p> [F + h}) * B . k * (F + h)* k - {p>F) * B hk * F* k ) e - p3( *- r) dr 



< 



M (l + |p| 2 ) {\W\\u + ||fc||„) {k\\FU\\B 3 . k \ * |fe|||„ 

+ \\B jtk *{F + h)\\ v \\h\\ v } (34) 

PROOF. The proof is similar to that of Lemma [l3[ except Lemma |l6| is used instead of Lemma O. 
□ 

Lemma 19 For F S a^cf v > Ap + 1 large enough so that (v/2)~P\\F\\ v < 1, AfF 
satisfies the following bounds 



\WF\\u < A b {T) ^Q(^)T( 3 ^)/ 3 (z,/2)-^ 



11*11 



W2)^||F||, 



(35) 



Further, {or h e such that (v/2)-P(\\F\\ v + < 1 

\\AT(F + h) - ATF\\ V 



\\h\\u 



< A b (T) VCi(0)(i//2)-°yT< a -'>/ 3 

" ^ ' [l-(«//2)-/»(||F||„- 



Proof. On inspection of ( |16[) and using Lemma |13|, it follows that 

3 oo 

ll^ll, < jc,.(0)T^)/ 8 E 111^1*1^111, P^ + ll^olL 



(36) 



(37) 



3=0 



k=0 



Now using Corollary |10| and noting the dependence of C\ on j through ctj, (35) follows. As far 
as (|36"|), from inspection of ( pif ) and using Lemmas [l?] and H, wc get 

+ - < ^CjiflTP-wU \B jfi \*\h\ \\ v 



3=0 



+ E (11*11" + WW WlWPi,*! * Nil- + * ( F + h)U\h\\„} (38) 



fc=i 
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On using Corollary 



10 



and noting the dependence of C\ on j through ctj, (|36|) follows. □ 



Comment 8: Lemma |T9| is the key to showing the existence and uniqueness of solution in 
to ([l6|), since it provides the conditions for the nonlinear operator M to map a ball into itself as 
well the necessary contractivity condition. 



Lemma 20 If there exists some b > 1 so that 
(y/2)-Pb\\F \\ v < 1 



3 



C J (0)(^/2)-^T( 3 -^/ 3 ^ 1 1 

3=0 

then the nonlinear mapping Af maps a ball of radius b\\Fo\\ u back into itself. Further, if 



^ rr^ fol <*-* ^ 



MT) U FPPW < 15 (40) 

then Af is a contraction there. 

Proof. This is a simple application of Lemma pi if we simply note that \\F\\* < b k \\Fo\\*. □ 



Lemma 21 For any given T > and cj> in the interval (0, 7r/6), for all sufficiently large v, 
there exists a unique F 6 A<f> that satisfies the integral equation ffl). 



Proof. 

We choose b — 2. It is clear from the bounds on ||fb||„ in Lemma 7 that for given T, since 
Oi r > 1) 6(^/2)~^||Fo|| ly < 1 for all sufficiently large v. Further, it is clear on inspection that 
both conditions ( p9|) and (jit]) arc satisfied for all sufficiently large v. The lemma now follows 
from the fixed point theorem. 
□ 

4.1 Behavior of solution F s near p = 

Proposition 22 For some K\ > and small p we have \F S \ < Ki\p\ Ur ~ 1 and thus \f s \ < 
K2\y\~ ar for some K2 > as \y\ — > 00 in T> p t p. 

Proof. Convergence in || • ||„ implies uniform convergence on compact subsets of K, and we can 
interchange summation and integration in (|l6|). With F s the unique solution of ( |l6| ) we let 

00 

G 3 = E(-i) 3 's ilfc * F: k 

and define the linear operator Q by 
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Clearly F s also satisfies the linear equation 

F s = QF S + F Q or F. = (1 - g)~ l F 

For a > small enough define S a — S H {p : \p\ < a}. Since F s is continuous in S we have 
lim a |o \\G\\ = 0, where the norm is taken over C(S a ). 

% @, (0)> © and Lemma |, we see that ||-F ||oo < # 3 |a| Qr ~ 1 in S a for some K 3 > 
independent of a. Then, as a j 0, we have 

max|F(p,t)| = ||F|| < (1 - HSU)" 1 max ||F || < 2^ 3 |a| a '- 1 

5 a S a 

and thus for small p we have i)| < 2i^3|p| Qr_1 and the proposition follows. 
□ 

Proof of Theorem ^. Lemma || shows that if / is any solution of (|l|) satisfying Condition 1 , 
then 6 ^0-5 for < <5 < for v sufficiently large. For large y, the series (^) converges 

uniformly and thus F = £~ 1 {f} satisfies (|l6|), which by Lemma |2~l| ha s a unique solution in 
A4, for any (ft G (0, 7r/6). Conversely, if € is the solution of (|l6|) for v > v\, then from 
comment 6, f s — LF S is analytic in V^ p for < (f> < <fi < tt/6, for sufficiently large p, where 
in addition from Proposition^, / s = 0(y~ ar ). This implies that the series in (Q) converges 
uniformly and by properties of Laplace transform, f s solves (1) and satisfies condition (1). 

Comment 9: Theorem |^can be applied directly to each of the examples 1-3 in the following 
sections to give existence and uniqueness of solution in £>0 lP for any given time T, provided p is 
large enough. However, this general theorem does not provide the specific dependence of p on 
T. In the following sections, we not only show that Theorem [| can be applied to the examples 
given, but use the specific information on bj,k(y,t) and r(y,t) to obtain dependence of p on T. 
This requires additional case-specific lemmas and theorems. 



5 Example 1 

This example comes in the context of solving the leading order inner-equation for a complex 
singularity of the conformal mapping function corresponding to an evolving Hele-Shaw flow 
(Jl2]|, Equations 5.5-5.9) as well as a 2-D dendrite in the small Peclet number limit, when surface 
energy anisotropy is small (]l7[, Equation A44 after a transformation). Consider the PDE 

Ht = H 3 H xxx (41) 

with the initial condition 

H{x,Q)=x 1 (42) 
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where < 7 < 1 (Note: this 7 is related to the (3 defined in , through 7 = /3/2), and the far 
field matching condition 



The transformations: 



H(x,t) = x"> + 0(x 4 ~>- 3 ) 

t 2-rr ' 

"3(1-7) ' 3(1-7)- 



(43) 



as |a;| — > 00 for arg x S (— 3( -?^ ■, , ) ). Here 7 is real and in the interval (0, 1). 



(44) 



1-7 

bring the equation to the form (Q), with 

r(y,t) = - 1 { 1 -l)- 2 ( 1 -2)y- 2 (45) 

and the only nonzero coefficients are {hi j : i, j = 0, 1, 2, 3}, are presented below as a matrix. 
227 - 11 7 2 - 



y z (7 - i) z 



n 67 - 37 2 - 2 50 7 - 25 7 2 - 18 (l- 2 7 )(2 7 -3) 



y i (7-ir 



2/ 5 (7 - 1) 



7 7 2 _ 14 7 + 6 77 2 - 14 7 + 6 „ 77 2 - 147 + ( 



(t-i)V 

-3JT 1 



y 3 (7 - i) 

-9y~ 2 



Sy" 1 
The initial condition (E2I) translates as 



y 4 (7 -i) 2 

-9y~ 3 

3y" 2 



2/ 6 (7 - 1) 
7 7 2 _ 14 7 _|_ 1 
y 5 (7- I) 2 

-3 3/ 



Also, the condition ( (43|) implies that as y — > 00 for (y, t) € 1^<p,p (as defined earlier), 

/(y,t) = 0(y- 2 ) 



(46) 



(47) 



From the expressions for bj t k above it is possible to calculate Bj^ explicitly. For our purpose, it 
is enough to note that aside from -63,0, which is identically zero, we can write 

\B hk {p,t)\<C\p\ 2 -J +k 

for a constant C independent of j and k, as well as T and therefore from Lemma 0, we conclude 
that for F € A4, for (j, k) ^ (3, 0), 



|||B^|*|F||U < C v- 3+j - k \\F\\ v 
We also note that from Lemma|7| that 

\\F (p,t)\\ u < A r v- x T 



(48) 



(49) 
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where A r here is independent of T. Since only a finite number of Bj t k are nonzero, it is better to 
use properties ( |37| ) and ( |3"8| ) directly to obtain conditions for the fixed point theorem to apply: 

3 3 

l + C J2'j2 bk A r T k ^ 2k +^ 3 T^^ 3 < 1 (50) 

j=0 k=0 



3 3 

Y^Y. bk ^ + A r T k iy- 2k+ ^ 3 T^/ 3 < 1 (51) 

Here the primes in the summation symbol in (|50|) and ( pl| ) mean that the term j = 3, k = 
is missing. Each of these conditions (|50| ) and (|5l|) are satisfied for j/ _1 T 1 / 3 sufficiently small 
for any choice of b > 1. The condition that T/v 3 is less than some number translates into 
T/p 3 being sufficiently small, i.e. we are restricted to a region of space in the x-plane where 
Ty x 3(i- T ) i s sufficiently small. It was noted earlier [Kadanoff, private communication, 1991 and 
independently by Howison Jl6| for a special case] that there was a similarity solution to ( pLT| ) of 
the form: 

H(x, t) = fr/Pfi-T) h (x/t^T)^j (52) 

The resulting ordinary differential equation for 71(77) was solved numerically Jl|] and a first few 
singularities (in order of distance from the origin) of H were thus determined. It was surmised 
that these solutions form an infinite set that straddles the boundary of the sector | argx| < 3?T3^ ; 
over which one can specify the asymptotic condition h(rj) ~ ?7 7 that one needs to satisfy initial 



and far-field conditions (42-43|). Later, these conclusions were confirmed rigorously by Fokas 
and Tanveer [l4| , who transformed the equation into Painleve Pn and used isomonodromic 
approaches for integrable systems to confirm the earlier behavior seen numerically. These results 
are also amenable to exponential asymptotic methods and formal trans-series association with 
actual functions, that have recently been worked out ||. The latter method is more general 
than the isomonodromic methods since they apply equally well for integrable and nonintegrable 
equations. The application of our existence and uniqueness results mean that the only solution 
to the initial value problem for the PDE are those with similarity structure given by (^2|) . 

From what has been discussed so far and proved in this paper, an interesting aspect of the 
complex plane initial value problem is that that the initial condition (^) is not recovered as 
t — > + , except in the sector | argx| < 3 /j_ 7 ) ■ This follows from the equivalence of small t with 



large x in the similarity solution (B2| 



6 Example 2 

The second example also comes from Hele-Shaw flow (|L2]], equations (6.10)-(6.12)) as well as 
dendritic crystal growth for weak undercooling and for weak anisotropic surface energy fl^| . In 
the asymptotic limit of surface tension tending to zero, it was determined that in the neighbor- 
hood of an initial zero, the local governing equation is 

H t + H x = H 3 H xxx — -H 3 (53) 
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The initial condition is: 

H(x,0) = x^ 1 ' 2 (54) 

The far-held matching condition with the "outer" asymptotic solution is 

H(x,t) = x- 1 ' 2 + 0{x- 5 ) (55) 

as \x\ — > oo with arg x 6 (— 4tt/9, 4tt/9). It is to be noted that in this case, unlike case I, there 
are no similarity solutions which satisfy both the initial and asymptotic boundary conditions. 
We introduce the transformation: 

x = t + (jfj ^ , Hix.^^x- 1 ' 2 + x-^ 2 y- 1 f(y,t) (56) 

Note that if x is large enough, y ~ | x 3 / 2 . The initial condition (|54|) implies that 

h(y) = (57) 
and the asymptotic far-field condition ( |55| ) implies: 

f(y,t) = 0(y- 4 / 3 ) (58) 
as y — ► oo in some 2?0, p . Under the change of variables, the PDE is of the form ([!]) with 

KM) = (59) 

and each 6 J; fe containing one or more terms of the form x~@ y -<5 , with /3 > and |/3 + 5 > 0. 
The exact expressions for 6^. are given in the appendix. 

Comment 10: Since the conditions for Theorem | hold, we may simply apply it and obtain a 
unique analytic solution over a sector T>^^ p for a fixed (f> satisfying conditions < 4> < 7r/6. 
The theorem, when applied to (|5^), would imply that for any T, there exists unique analytic 
solution H(x,t) in the sector arg x € (— 7r/3— 30/2, — tt/3 + 3</>/2) provided |x| is large 
enough. Since this is true for any (j> in the interval (0, 7r/6), the theorem establishes the existence 
of a unique analytic solution assumed before fl2f . However, the restriction on how large x has to 
be depends on T and because of the generality, Theorem || does not give a precise dependence 
on T. Finding this constraint is the objective of the rest of this section. 

Lemma 23 Let g{y,t) — x~P y~ 5 , where x is given by (f5^), (3 > and |/3 + <5 > 0, then 
for any p <E S^, 

\G(p,t)\ < C piP+z- 1 (60) 

where C is independent ofv and T, but can depend on (3 and 6. Also, if in addition, |/3 + <5 > 1, 
then for v > 1, 

||G||„ < Cv-^l 3 -^ 1 (61) 
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Proof. Given the relation between x and y, we note that we may write 

g{y,t) = y-^- s h(ty-^) 

where 

g\ -2/3/3 / 2 2/ 3 \ -0 



It is to be noted that for arg s bounded away from ±7r, h(s) is uniformly bounded. It is also 
clear that 

G(p,t) = C^gfat] = ^~P 2P/3+S ^ (J e s h{t P 2 ^s- 2 ^) ds\ 

where the contour C is similar to that shown in Fig. 1, except in the s-plane. The intersection 
point of C on the real s-axis will be chosen to be 1. It is clear that for p G S$ and s on the 
contour C, arg(ip 2 / 3 s~ 2 / 3 ) S (— 5n/9, 5n/9). Therefore h is uniformly bounded and therefore 

/>oo 

\G{p,t)\ < CH^"- 1 / e - r ' 2 dr 



Hence the lemma follows. □ 

Comment 11: Note that the preceding Lemma gives a much sharper result than applying the 
more general Lemma 3, because we made specific use of the form of the function g(y,t). 

Corollary 24 \\R\\ V < C v- 1 / 3 . 

Proof. This follows simply from the expression for r and and application of Lemma □ 
Corollary 25 For some C independent of T and p, 

(i) |Bo.o| < C \ P \ 2 , |S ,il < C H 4 / 3 , |B 0l9 | < C \p\ 3 and \B , 3 \ < C \p\ 7 

(ii) |fli,o| < C\p\, < CM 8 / 3 , \B h2 \ < C\p\ 13 / 3 and \B 1>3 \ < C\p\ e . 
(w)\B 2 ,o\ < C,\B 2 ,i\ < C \p\ r ° /3 , \B 2 .2\ < C H 10 / 3 and \B 2 . 3 \ < C \p\ 5 . 

(iv) \B S>0 < CT |p|- 1/3 , \B 3ll \ < C \p\ 2 / 3 , \B 3 , 2 \ < C \p\ 7 / 3 and \B 3 , 3 \ < C \p\ 4 . 

Proof. These inequalities follow immediately from the expressions of the coefficients bj t k in the 
appendix and Lemma E3|. □ 



Corollary 26 For C independent of T ', p, j and k and v > 1, 

1115^1*1^11, < C ^79-2/3-^/3 ll^l^ for(i)fe) ^ (3)0) 

IH^olH^IIU < CTv- 2 ! 3 \\F\\ U 

Proof. The estimates follow immediately on examination of the upper-bounds on B 3 k and 
using Lemma ^ with p = 0. □ 

Comment 12: Application of Lemma || to the results Corollary leads to stronger bounds; 
however the bounds indicated in Corollary Efj] suffice for our purposes. 
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Lemma 27 For any < (j> < 7r/6, the integral equation (|7qj, urafft and R as determined 
in this section, has a unique solution F s G A^ provided T / 3 < e, where e ( depending only 
on <j>) is small. Further, F s = 0(p 1 ^ 3 ) as p — > in 

PROOF. Since F„ (p, t) is given by @ , and Fj (y) = 0, it follows that ||F 1 | v < CTz/" 1 / 3 . 
Further, from the estimates of Corollary ^6|, it follows from (^) that the condition for mapping 
a ball of radius &||i*o||i/ back into itself is now given by 

\ + Cj2' (Tv 2 ' ' (bTv-^Y + CTv-^ < 1 

j,k=0 

where J^' denotes is the summation without (j,k) = (3,0) term. From (|38|), the contractivity 
requirement becomes 

3 

C ( k + bk (Tv~ 2/3 ) ° (i>IV~ 2/3 ) + CT^ 2/3 < 1 

j,k=0 

It is clear that both conditions are satisfied for some b > 1 (say b = 2) if T V - 2 ^ is chosen 
sufficiently small. This ensures existence of a solution in the Banach space A^. Since C in the 
above equations depend on (f>, the upper bound of T v~ 2 ^ for which the solution is guaranteed 
to exist depends on <j>. Further, since r(y, t) = 0(y~ 4 / 3 ), by applying proposition |22^, it follows 
that F s (p,t) = 0(p 1 / 3 )a,sp -> in S . 
□ 



Theorem 28 For any < < 7r/6, i/iere exists a unique solution H(x,t) satisfying (pqj-fp, 
that is analytic in x in the domain 

( 4 2 4 2 \ 
(x,t) : |x| > p , arg a: G I --tt + -</», -vr - -0 J , <t<T 

provided TpT 1 < e, with e (depending only on <f>) small enough. Thus, for any t, when 
arg x G (— 47r/9, 4tt/9), there exists a unique analytic solution satisfying ( jj^p^ j when t/\x\ is 
small enough. 

Proof. By applying the equivalence between the solution to the integral equation (fig), F s G A<j, 
for < <f> < 7r/6, to the analytic solution f s of the partial differential equation (|IJ) satisfying 
condition 1 (as shown in the proof of theorem 2), it follows that in this particular example, a 
solution f s exists for (y, t) G provided T p~ 2 / 3 < e and that f s = 0(y -4 / 3 ) as y — » oo 

in 2?0 iP . Theorem 28 follows merely from noting the change of variable (y, t, /) to (x, t, H), once 
we choose p — p^ 3 and use the relation y ~ |x 3 / 2 for large x. □ 



7 Example 3: Strongly anisotropic inner equation 

For strong anisotropic surface energy, the analytically continued conformal mapping function 
that maps the upper half-plane to the exterior of a one-sided two dimensional dendritic interface 
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for small Pcclct number satisfies, upon transformation, the following leading order inner-equation 
near a singularity of particular type (See equation (A16) after some elementary transfor- 
mations): 

H t = ff 1/3 H xxx (62) 

with initial and boundary conditions 

H(x,0) = x~ 95 (63) 

H(x,t) = x- gs + 0(x' 125 ~ 3 ) as \x\ -> oo (64) 
for arg x € ^ — 3(TT57 ' 3(i+<5) ) wnere <5 > 0. If we introduce the transformation 

y = YTs i = x ~* s ( l + fMlv) (65) 



then we obtain an equation of the form ([!]) with 



, , 95(95 + 1) (95 + 2) 
rM = \, + 1 )V (66) 



, 9J(9J + 1)(9^ + 2){(1 + //y) 4 / 3 - l}t/ 

bo{f,y,t) = rrr^m 

{5 + 1)V/ 

+ (5A6 2 + 277S + 32) y J J ' y J _ (67) 
(0 + l) z jr 

" (-W-^- 6 ) (1 + //rf,V2 (68 » 

w . r , t) . ^±ifct(Z^ (69) 

2/(<J + 1) 

6 3 (/,y,t) = 1 - (l + f/y) 1/3 (70) 
Using the series expansions in f/y it follows that 

95(95 + 1)(95 + 2)(# 3 1 ) (545 2 + 2775 + 32)^) 



6, 



'(),*■ 



(5 + l) 3 l/ 4+& ((5 + l)2y3+fc 



. 2175 + 26 48,5 1 A/3\ _ k _ 2 



6 



2. A: 



3(95+l)( 1/3 



k 



y k+1 (6 + l) 
h, = 

b 3 . k = -f^ 3 ) forfc ^ 1 
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It is clear that 

Bo,k 



95(95+l)(9S + 2)(^\) P k + :i (54 ( ? 2 + 277J+32)( 1 /> fc + 2 
(<5 + l) 3 (fc + 3)! + (<5 + l) 2 (fc + 2)! 

217(5 + 26 48(5 1 /l/3\ p k+1 



B 



2./v 



(l + <5) 2 (5 + 1 J V fc / (fc + 1)! 
3(9(5+ 1)( 1 / 3 )/ 



((5 + l)fc! 
B3.0 = 

Thus, the following estimates hold: 

111^*1*1*111,, < Cv- k +t-* \\F\\ V 

where C is a constant that can be made independent of j, k, T. Therefore, using the above 
relation, from (|37]), the condition for mapping a ball of radius b \\Fq\\ u back into itself becomes: 

3 00 1 
C^'^ v ^ a ~ i)/S ( v ~ lb \\Foh) h + \ < 1 

j=0 fc=o 

where the ^2' indicates that j = 3, k = term is missing from the summation. Applying the 
estimates of this section to (p8[), the contraction condition is: 

3 00 

^E'E(^ 3T ) (3 ^ )/3 ( fc + 1 )K 1& n F o^) fe < 1 

Since ||i*b||i/ < KTv^ 1 , a sufficient condition for use of contraction mapping theorem is that 

bKTv- 2 < 1 

i.e. that Tv~ 2 is small enough. Note that in that case Tv~ z is automatically small when v 
is sufficiently large. The restriction Tv~ 2 small means that the differential equation (|62|), with 

conditions ( |63[ - p4| ) has a unique analytic solution for any x 6 ( — 3(1^5) ' 3(i+<?) ) * n a re §i° n 

where when tx~ 2S ~ 2 is small enough. However, ( |62] ) admits a similarity solution 

and (7(77) solves an ordinary differential equation and the asymptotic boundary condition 

2tt 2tt 



q(rj) ^ rj 96 for 7] — * 00, arg 77 6 



"3(1 + <5)' 3(1 + 5) 



Uniqueness means that the similarity solution is the only solution. However, the restriction 
££-2(5-2 <^ x i s sub-optimal; from the similarity structure of the solution, one expects analyticity 
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in a sector for tx~ 35 ~ 3 <C 1, i.e. the integral equation (16) should have a unique solution for 
TV~ 3 sufficiently small. The condition for convergence of the infinite series involving the norm 
of F is stronger than the original condition for the convergence of the infinite series in / that 
appear in the expansion of (1 + //y) 1 ^ 3 . This is the reason for the sub-optimal estimates in this 
example involving infinite series. Nonetheless, the uniqueness results even for a restricted range, 
prove that the similarity solution determined earlier is the only solution to the problem. 



8 Conclusion 

We have proved existence and uniqueness of solution to a class of third-order nonlinear partial 
differential equation in a sector of the complex spatial variable y for sufficiently large Our 
technique is akin to Borel summation, that has demonstrated its effectiveness in the analysis of 
general classes of nonlinear ODEs ||. The class of PDEs for which existence and uniqueness 
has been proved contains three examples that arise in the context of Hele-Shaw fluid flow and 
dendritic crystal growth. The uniqueness results shows that the similarity solutions assumed 
earlier for examples 1 and 3 are the only ones that satisfy the given initial and far-field matching 
conditions. Accordingly, the singularities of the PDE solutions are those that correspond to 
singularities of the similarity solutions. 
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Expressions for bj t k for Example 2 

bo,0 = —— o 7" ^ 

6(),1 = — 
6o,2 = — 



bo,s 



62,1 
62,2 
62,3 



63,2 
63,3 



35 1 75 ^_ 9/2 _ 45 (12) 1/3 _ 15 (18) 1/3 



6 x 3 / 2 y 2 4 8 x 7 / 2 y 2 / 3 4 x 5 / 2 y 4 / 3 

35 1 45 3 45 (18) 1/3 135 (12) 1/3 



2 x 5 / 2 y 3 x 11 / 2 y 2x 2 y 4 x 7 ' 2 y 7 / 3 8 x 9 / 2 y 5 / 3 
165 1 135 (12) 1/3 1 35 1 45 (18) 1/3 



4 x 13 / 2 y 2 8 x 11 / 2 y s / 3 2x 3 y 2 2 x 7 / 2 y 4 4 x 9 / 2 y 10 / 3 

45 (12) 1/3 35 1 105 1 15 (18) 1/3 

Y 3:13/2^11/3 _ Y a^/V 8~ i«/y ~ T a; 11 /2 2y i3/3 



15 (18) 1/3 45 ^y(12) 1/3 35 1 
1,0 " 4 a; 5 /2^ + 8 x 7 1 2 + 6 :r3/ 2 y 

35 1 45 (18) 1/3 135 (12) 1/3 
61,1 — r. /■> ■> + —r 



2 x 5 / 2 y 2 4 i 7 / 2 ?/ 4 / 3 g x 9 / 2 y 2 / 3 

35 1 135 (12) 1/3 45 (18) 1/3 
01,2 = ; H ~ n ,„ „>, + —r 



2 a; 7 / 2 ?/ 3 8 a: 11 / 2 ^ 4 ^9/2^7/3 
35 1 15 (18) 1/3 45 (12) 1/3 

01,3 = — Km * + ~. 11 ;o i.« + TT 



6 a; 9 / 2 ?/ 4 4 a^/Zyio/s 8 x 13 / 2 y 8 / 3 
,2,0 = -3^-^ 



4a .5/2 

9 _ 27 (18) 1/3 
' x 5 ' 2 y ~ T a- 7 / 2 ^7/ 

9 27(18) 1/3 
" x 7 / 2 y 2 ~ Ax 9 / 2 y 4 / 3 

3 9(18) 1/3 



63,0 = -1 + 



x 9/2 y 3 4a.ll/2y7/3 

3y 



1 9 _ 5 /2 

63,1 = -Z 



2:r 3/2 

2 
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2a- 7 / 2 ?/ 
3 

2x 9/ 2 y 2 
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